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Waveguide Fiber 
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(Manuscript received September 1, 1982) 

In this paper we compare modal analyses of a Modified Chemical Vapor 
Deposition (MCVD) optical waveguide fiber and an idealized fiber. The 
pertinent profile parameters of the idealized fiber are obtained from the MCVD 
profile. We have validated the concept of an optimum a developed in our 
previous work. The reduction of the bandwidth of the MCVD fiber is shown 
to be directly related to the imperfections in the profile inherent in the MCVD 
manufacturing process. 

I. INTRODUCTION 

Multimode lightguides, produced by the Modified Chemical Vapor 
Deposition (MCVD) process, are expected to be of continuing impor- 
tance to the Bell System. Such lightguides with bandwidths exceeding 
1 GHz • km have been produced in the manufacturing plant. 

The unique MCVD process produces the index gradient in the core 
by depositing many layers of material with different compositions, 
leading to a ripple effect in the profile. Other imperfections in the 
MCVD process include certain mechanical limitations, difficulties at 
the core-cladding interface, and the altering of the core center index 
(burn-off) during the collapse of the preform. These imperfections in 
the MCVD process lead to a unique modal structure in the fiber, which 
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may differ radically from the modal structure obtained from a perfect 
power law profile. These differences, for a particular fiber preform and 
idealization, are analyzed and compared. The effects of the imperfec- 
tions on dispersion and bandwidth are evaluated. 

II. EXPERIMENT 

Our method 1 for solving Maxwell's equations has been described 
earlier and we shall not repeat it here. A brief schematic of the process 
is given in Fig. 1. In the experiment, the perfect power law profile (N{) 




Fig. 1 — N-, is a perfect power law profile. N* is a typical MCVD profile that must be 
normalized in a particular way to be compatible with the calculations described in Ref. 
1. From these calculations we obtain the G function, from which one further determines 
the effective indices, N e , and the group indices, N g . Further processing leads to the 
calculation of the Poynting vector, the field functions, the delay times, etc. 
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is replaced with the actual MCVD profile (N x ), which differs from the 
perfect profile. Thereafter, the data are processed identically to find 
the effective indices N e , group indices N e , the Poynting vector, and 
other propagation parameters. 

An MCVD preform was given to us by S. Jang of the Western 
Electric Research Center in Princeton, N.J. This preform had been 
used earlier to produce fiber for experiments. It was typical in the 
sense that the bandwidths averaged 1 GHz -km. By means of the laser 
beam refraction method 2 the refractive index profile was determined 
and is shown in Fig. 2. 

We note in this profile the familiar characteristics of the MCVD 
preform, i.e., a burn-off region in the center, an undefined core- 
cladding interface region, and ripples on the profile, which are more 
pronounced as we near the core center. 

The pertinent data required for our analytical testing are the ap- 
proximate a and the maximum AN. These data can be obtained from 
known procedures, 2 or alternatively, as in our case, by utilization of 
the technique described in one of our earlier publications. 3 Briefly, at 
a selected wavelength, one obtains the modal display for zero azimuthal 
numbers, fits the data to a linear least squares equation for iV g vs. iV e , 











0.016 


"7~^A 






0.012 


- 














^ 0.008 








0.004 


- 




\ / 








C--JI 


-0.004 




I I 


I 



0.08 0.16 0.24 

RADIUS IN CENTIMETERS 



Fig. 2— Typical MCVD profile before normalization. 
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and compares this slope to the slope vs. a, which has been previously 
calculated for that wavelength. The parameters to be used in this case 
were a = 2, AiVis typical for an 11.5-percent Ge(>2 dopant, and radius 
{R cc ) = 25 urn. 

These data are then used to produce modal displays for perfect 
power law profiles (IDEAL) at the wavelengths and polar indices used 
for probing the experimental (MCVD) profile. The programs are run 
in parallel, to be used eventually for comparison purposes only, i.e., 
there is no communication between the two series of calculations. 

Inherent in our calculations is the automatic detection of significant 
mode splitting between the EH vs. HE, or TE vs. TM, and in the 
displays to follow, the splittings are shown, wherever possible. 
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Fig. 3 — Composite data for X = 0.9 fim. The top portion (a) is the IDEAL. The 
bottom portion (b) is the MCVD. The straight lines are the linear least-squares curves 
fitted to all the data in each case, exclusive of those modes near the cladding, and for 
those modes widely divergent from the line. In general, the open circles are for polar 
indices m = 1, 3, 5, • • • last, while the solid circles are for polar indices m = 0, 2, 4, 
last. Split modes are not readily identifiable in this display. 
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Fig. 4— Same as Fig. 3, A = 1.0 urn. 
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Fig. 5— Same as Fig. 3, X = 1.1 ^m. 
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Fig. 6— Same as Fig. 3, X = 1.2 urn. 



We selected eight wavelengths (X) spanning the region 
0.9 < X < 1.55 mhi. 

Each wavelength was checked for all possible modes (m = 0, 1, 2, 
• • • last) for both the IDEAL and the MCVD profiles. Because the 
extremely large amounts of data required in this analysis involved 
exorbitantly large amounts of computer time, we reduced the required 
accuracy to a minimum sufficient for our purposes, thereby saving 
considerable time and expense. Nevertheless, much computational 
time was required to accumulate the data given in this report. 

Our results on the foregoing procedure are displayed in Figs. 3 
through 10. In each figure, the dots represent the modes actually 
determined, the line is the linear least-squares curve fitted to the data, 
excluding some of the modes near cutoff, the upper portion (a) is for 
the IDEAL profile, and the lower portion (b) is for the MCVD profile. 
Again, some modes that are well off the chart are not shown because 
of space limitations and were not used to calculate the fitted curves in 
any case. All the fitted curves, translated where necessary but not 
rotated, have been condensed into the single plot given in Fig. 11. 
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III. DISCUSSION 

The degree of scatter in the modes for the IDEAL case (Figs. 3 
through 10, Part a) is due in part to the inherent splitting of the 
HE-EH modes, which has been further compounded by the reduced 
accuracy used in this analysis. 

The scattering of the modes for the MCVD case (Figs. 3 through 
10, Part b) is much greater than the IDEAL series, perhaps by a factor 
of 5 to 10. Obviously, the greatly increased scattering for the MCVD 
series is due to the imperfections in the MCVD profile, such as burn- 
off and ripple. The contribution of each of these imperfections to 
increased modal dispersion (poor bandwidth) has been determined 
and will be the topic of a forthcoming paper. 

Further, within the statistical precision and accuracy employed by 
the authors, one can calculate that the IDEAL and MCVD data will 
both extrapolate to the same index at the core. We should realize that 
the MCVD index of refraction data of Fig. 2 was obtained at X = 
0.6328 /xm. This must be converted to the proper AN at the probing 
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Fig. 8— Same as Fig. 3, A = 1.32 /im. 



X's used in any analysis, and was done in this analysis. The results 
indicate that our parameter AN as calculated at each X is reliable. 

The parallelism of the least-squares fitted lines, seen in Fig. 11, 
indicates that our calculations of the parameter a from the MCVD 
profile are also reliable, and that an idealized profile may be used to 
predict characteristics of an MCVD profile. We must realize that a 
does not vary with X, and this is true regardless of whether we are 
using the IDEAL or MCVD profile. We are immediately struck by the 
rotation of these curves as we change X. We have noted a similar 
effect 3 at a fixed X for a changing a. It should be noted that the 
previous work 3 and the present work agree closely. 

We have one further observation in the present work, which is not 
readily seen in Figs. 3 through 10: although both the IDEAL and 
MCVD profiles were normalized to a 25-/tm radius, the IDEAL profile 
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Fig. 9— Same as Fig. 3, X = 1.40 nm. 



generated more modes than the MCVD, without exception. To account 
for the fewer MCVD modes a crude estimation of the effective radius 
for the MCVD indicates that 

r t . ~ 23 Mm. 

It may be only coincidence that this difference of 2 /*m is of 
approximately the same magnitude as the width of the burnout. 

IV. REDUCTION OF DATA 

Two useful measures of dispersion were derived for this report. The 
first, which is relatively simple, is the calculation of the dispersion due 
to the rotation of the line as a function of A. This implies that every 
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Fig. 10— Same as Fig. 3, X = 1.55 fim. 

group index (N e ) would be on the line. However, the total spread in 
N e is now a function of the slope of the line, i.e., the spread in N e 
between the HE ia and the last possible mode that can be obtained on 
that line. We have given this measure the name "sigma of rotation", 
(or). Of course, we must count the number of modes at each N e to 
find the mean N g and then calculate or in the conventional manner. 

The second measure of dispersion we have utilized is the standard 
deviation of the group indices determined from the indicated least - 
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Fig. 11— Composite data. Least-squares fitted lines (IDEAL and MCVD) from the 
data of Figs. 3 through 10, identifiable in this display. 



squares fitted line for N g vs. N e . Thus, if the curve value of N e is 
denoted by JV K (curve), then 



V 2[N E - iV g (curve)] 5 



where n is the total number of modes. 

The results of these calculations for o- R are given in Fig. 12, and as 
indicated in Fig. 11, are about the same for the IDEAL and MCVD 
fibers. Further, if this were the only measure of dispersion employed, 
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Fig. 12— o R vs. X for IDEAL or MCVD. 



we would expect this fiber to yield the maximum bandwidth at A = 
1.00 /im. This agrees precisely with our previous work. 3 The results 
for the calculation of a Ng are given in Fig. 13 and we can see that: 

o„„(MCVD) * 10 o».(IDEAL). 

Also, we note further that o Ni is a decreasing function for larger values 
of X, and becomes virtually constant when X > 1.2 itm. 

A better understanding of the roles of the two measures of dispersion 
previously described is given in Figs. 14 and 15. In Fig. 14, or and o/v g 
are plotted for the IDEAL case, and in Fig. 15 the same data are 
plotted for the MCVD case. Apparently, the o that is dominant at any 
X will be the principal cause of inability to maximize the bandwidth. 

From theory 4 we have calculated the total dispersion o c by: 

o c = v Or + On v 

and these data are given in Fig. 16. We see that in the region where 
or is dominant (X > 1.2 //m), the curves are roughly parallel. The 
larger o c for the MCVD in this region comes about because the a Ng of 
the MCVD is greater than a Ng of the IDEAL. 
Over the range 0.9 /zm < X < 1.2 /mi, the IDEAL case peaks (o- c = 
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Fig. 13— <T/v, vs. X for IDEAL and MCVD. 
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Fig. 17— l/2ira c for the IDEAL case. 

minimum) sharply at about 1.0 /tm, while the MCVD displays a very 
broad response that appears to peak at about 1.1 pan. A qualitative 
measure to estimate the pulse shape can be realized by calculating 
W(h/2); W(h/2) = 1/2tt<x c (MHz). 

These data are given in Figs. 17 and 18 for the IDEAL and the 
MCVD, respectively. The IDEAL gives a bandwidth of 9 to 10 GHz/ 
km at 1.0 /im and falls off to half-height at ±0.7 urn. The MCVD 
obtains a maximum of about 1.5 GHz km at «1.1 /tm, and falls off to 
half height at ±0.2 ^m. 



V. CONCLUSION 

There can be no doubt that the concept of an optimum a is valid 
for perfect power law and for experimental (manufactured) profiles. 3 
The reduction of the bandwidth is directly related to the imperfections 
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Fig. 18— 1/27t«t c for the MCVD case. 

of the MCVD profile such as burnout, ripples, etc. The direct role of 
each of these deviations from an idealized profile, in the reduction of 
bandwidth, has been calculated, and will be the subject of a forthcom- 
ing paper. These observations have been substantiated in many re- 
spects from direct communication with the materials scientists at Bell 
Laboratories in Murray Hill and Atlanta, and the Western Electric 
Company Engineering Research Center. We have also ascertained that 
we can improve the calculation of <r Ng by about a factor of 2. Thus we 
find the bandwidth for the idealized fiber is about 18 GHz • km, while 
for the MCVD fiber it is about 2.8 GHz -km. These bandwidths can 
be further improved by careful computer-aided design of the profile, 
for which investigations are currently in progress. 
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